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Using the coupled cluster method for high orders of approximation and complementary exact 
diagonalization studies we investigate the ground state properties of the spin-1/2 J\-Ji frustrated 
Heisenberg antiferromagnet on the square lattice. We have calculated the ground-state energy, 
the magnetic order parameter, the spin stiffness, and several generalized susceptibilities to probe 
magnetically disordered quantum valence-bond phases. We have found that the quantum critical 
points for both the Neel and collinear orders are J 2 cl w (0.44 ± 0.01) Ji and Jf w (0.59 ± 0.01) Ji 
respectively, which are in good agreement with the results obtained by other approximations. In 
contrast to the recent study by [Sirker et al. Phys. Rev. B 73, 184420 (2006)], our data do not 
provide evidence for the transition from the Neel to the valence-bond solid state to be first order. 
Moreover, our results are in favor of the deconfinement scenario for that phase transition. We also 
discuss the nature of the magnetically disordered quantum phase. 



I. INTRODUCTION 

Quantum phase transitions between scmiclassical mag- 
netically ordered phases and magnetically disordered 
quantum phases which are driven by frustration at- 
tract much interest; see, e.g., Ref. [l|. A canoni- 
cal model for studying such transitions is the spin- 
1/2 Heisenberg antiferromagnet with nearest-neighbor 
J\ and frustrating next-nearest-neighbor J 2 coupling 
{J\-Ji model) on the square lattice. This model 
has attracted a great deal of interest during the last 
20 years (see^ _R^s^J201ilIMM3IMMl3liIllI2IIl^ 



references therein). Recent interest in this model 
comes also from the synthesis of layered magnetic 
materials Li 2 VOSi0 4 , Li 2 VOGe0 4 , VOM0O4, and 
BaCdVO(P0 4 ) 2 (Refs^^i^i^) that might be described 
by the J1-J2 model. A new promising perspective is also 
opened by the recently discovered layered Fc-based su- 
perconducting materials^ 6 - which may have a magnetic 
phase that can be described by a J\-Ji model with spin 
quantum number s > l/2 i 37 ' 38 ' 39 

For the square-lattice spin-1/2 J1-J2 model it is well 
accepted that there are two magnetically long-range or- 
dered ground state (GS) phases at small and at large 
Ji separated by an intermediate quantum paramagnetic 
phase without magnetic long-range order (LRO) in the 
parameter region Jf 1 < J 2 < J 2 2 , where J^ 1 « OAJi and 
Jf 2 « O.6J1. The magnetic phase at low J 2 < J^ 1 ex- 
hibits semiclassical Neel LRO with a magnetic wave vec- 
tor Q = (tt, tt). The magnetic phase at large J 2 > J| 2 
shows so-called collinear LRO. It is twofold degener- 
ate and the corresponding magnetic wave vectors are 
Qi = (tt, 0) or Q 2 = (0,7r). These two collinear states 
are characterized by a parallel spin orientation of near- 
est neighbors in vertical (horizontal) direction and an 



antiparallel spin orientation of nearest neighbors in hor- 
izontal (vertical) direction. 

The nature of the transition between the Neel and 
the quantum paramagnetic phases as well as the prop- 
erties of the quantum paramagnetic phase and the 
precise values of the transition points, is still under 
debat e 2 - 3 - 4 - 5 - 6 - 7 - 8 - 9 - 10 - 11 - 12 ' 13 - 14 ' 15 - 16 ' 17 - 18 - 19 - 20 - 21 - 22 - 23 - 24 - 25 - 26 - 27 - 28 ' 29 - 30 ' 

In particular, stimulated by the recent discussion of 
deconfincd quantum criticality in two-dimensional 
spin systems ; 40 ' 41 a renewed interest in the nature 
of the phase transition between the semiclassical 
Neel phase and the quantum paramagnetic phase 
has emerged ! 27 ' 30 ' 42 ' 43 However, in spite of numerous 
intensive efforts focused on the transition between the 
Neel and the quantum paramagnetic phases in the 
J\-Ji square-lattice antiferromagnet and some other 
candidate models ; 44 ' 45 ' 46 ' 47 ' 48 ' 49 this field remains still 
highly controversial. 

For completeness we mention that the classical square- 
lattice J\-Ji model (s — > 00) exhibits a direct first- 
order transition between Neel state and collinear state 
at J 2 /Ji = 1/2. 

Recently, several extensions of J\-Ji model have been 
studied. Interestingly, with increase in the space di- 
mension from D = 2 to D = 3 the intermediate 
quantum paramagnetic phase disappears! 29 ' 50 ' 51 Also 
spatia l 52 ' 53 ' 5 ^ 5 -' 5 ^ and spin anisotropie a 26 ' 57 '^ 8 - as well as 
the spin quantum number s (Refs. [2l l30ll59l andlfjOT) have 
a great influence on the GS phase diagram. 

The goal of this paper is to study the GS phase di- 
agram for spin-half Ji-J 2 model on the square lattice 
using a high-order coupled cluster method (CCM). We 
complement the CCM treatment by exact diagonaliza- 
tion (ED) for finite lattices for a qualitative check of our 
CCM data. By calculating GS quantities such as the en- 
ergy, the magnetic order parameter, the spin stiffness and 



2 



generalized susceptibilities we will investigate the quan- 
tum phase transitions present in the model as well as 
the properties of the quantum paramagnetic phase. We 
will compare our results with the ones obtained recently 
using series expansions^ 

The CCM, introduced many years ago by Coester 
and Kiimmel^i is one of the most universal and most 
powerful methods of quantum many- body theory. For 
a review of the CCM see, e.g., Ref. |62|. Starting in 
1990 it has been applied to quantum spin systems with 
much SU ccess.^^^^^^^^^^^.^^^i^ A main 
advantage of this approach consists in its applicability 
to strongly frustrated quantum spin systems in any di- 
mension. With the implementation of parallclization in 
the CCM cod o 72 ' 76 high-order calculations are now pos- 
sible (see Sec.|ll|, improving significantly the accuracy in 
the investigation of quantum phase transitions driven by 
frustrationM^SLSmiiffi 

The Hamiltonian of the considered J1—J2 model reads 

h = j x SjSj + j 2 Sjs.,- , (1) 

iv) [ij] 

where J\ is the nearest-neighbor exchange coupling and 
J2 is the next-nearest-neighbor exchange coupling. Both 
couplings arc antiferromagnetic, J\ > and J2 > 0. In 
our CCM and ED calculations we set J\ = 1. We consider 
spin quantum number s = 1/2, i.e., s? = 3/4. 

The remainder of the paper is organized as follows. In 
Sec. |IT] we briefly discuss the CCM approach and illus- 
trate how to calculate of GS quantities of spin model (fTJ) . 
We present our results for the GS energy, the magnetic 
order parameter and the spin stiffness in Sec. IIII Al In 
Sec. lIII 51 we consider in more detail the phase transition 
between the Neel state and the quantum paramagnetic 
state as well as various susceptibilities testing the nature 
of the nonmagnetic phase. Finally, in Sec. IIVI we sum- 
marize our findings. 

II. COUPLED CLUSTER METHOD 

We start with a brief illustration of the main features 
of the CCM. For a general overview on the CCM the 
interested reader is referred, e.g., to Refs. 
and [Zllzizl. The starting point for a CCM calculation 
is the choice of a normalized model (or reference) state 
$), together with a set of mutually commuting multispin 
creation operators Cj which are defined over a complete 
set of many-body configurations /. The operators Cj are 
the multispin destruction operators and are defined to be 
the Hermitian adjoint of the C 1 /. We choose {|<I>);C/} 
in such a way that we have ($|C/ = = Cj\&), VI 7^ 
0. Note that the CCM formalism corresponds to the 
thermodynamic limit N — > 00. 

For the spin system considered, for |$) we choose the 
two-sublattice Neel state for small J2 but the collincar 
state for large J2. To treat each site equivalently we 



perform a rotation of the local axis of the spins such that 
all spins in the reference state align along the negative z 
axis. In the rotated coordinate frame then we have |$)= 
I |)| |)| I) ■ ■ ■ and the corresponding multispin creation 
operators then read Cf = sf , s^Sj", sfsjs^, where 
the indices k, . . . denote arbitrary lattice sites. 

The CCM paramctcrizations of the ket- and bra- 
ground states are given by 

H\V} = E\V}; (&\H = E(V\; 
|v&) = e s |<I>), S = £s,C+; 

IjtO 

(*| = ($|5e.- s , S = 1 + J2SiCT- ( 2 ) 

The correlation operators S and S contain the correla- 
tion coefficients Si and Si that have to be determined. 
Using the Schrodingcr equation, H\^f) = E\^f), we can 
now write the GS energy as E — (<f>|e" s 'i/e 5 |$). The 
magnetic order parameter is given by 

1 N 

M=--£(*| S ?|*), (3) 

i=l 

where s\ is expressed in the rotated coordinate system. 
To find the ket-state and bra-state correlation coefficients 
we require that the expectation value H = (\E , |.H"|\I') is 
a minimum with respect to the bra-state and ket-state 
correlation coefficients, such that the CCM ket- and bra- 
state equations are given by 

(<S>\Cye- s He s \<f>) = 0, V/ ^ 0, (4) 
($\Se- s [H, C+]e s |$) = 0, VI ^ 0. (5) 

Each ket- or bra-state [Eq. [4j or (5] belongs to a particu- 
lar index / corresponding to a certain set (configuration) 
of lattice sites k, . . . in the multispin creation opera- 
tor Cf = sf, sfsj~, sfs~jSfr, • • • ; see above. 

Though we start our CCM calculation with a ref- 
erence state corresponding to scmiclassical order, one 
can compute the GS energy also in parameter regions 
where semiclassical magnetic LRO is destroyed, and it is 
knowit^^^2JZJl that thc CCM yields p rec i S e re- 
sults for the GS energy beyond the transition from the 
semiclassical magnetic phase to the quantum paramag- 
netic phase. The necessary condition for the convergence 
of the CCM equations is a sufficient overlap between the 
reference state and the true ground state. 

It has been recently demonstrated 7 ^ that the CCM can 
also be used to calculate the spin stiffness p s with high ac- 
curacy. The stiffness measures the increase in the amount 
of energy when we twist the magnetic order parameter of 
a magnetically long-range ordered system along a given 
direction by a small angle 9 per unit length, i.e. 
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FIG. 1: Illustration of the twisted reference states used for the 
calculation of the spin stiffness p s . The angles at the lattice 
sites indicate the twist of the spins with respect to the Neel 
or the collinear state, (a): Twisted Neel state; the twist is 
introduced along rows in x direction, (b): Twisted collinear 
state; the twist is introduced along rows in e x + e y direction. 



where E(9) is the GS energy as a function of the imposed 
twist, and N is the number of sites. In the thermody- 
namic limit, a positive value of p s means that there is 
magnetic LRO in the system, while a value of zero reveals 
that there is no magnetic LRO. To calculate the spin stiff- 
ness within the CCM using Eq. ([6]) we must modify the 
corresponding reference states (Neel or collinear) by in- 
troducing an appropriate twist 6, see Fig. [TJ Thus the 
ket-state correlation coefficients Si [after solving CCM 
equation ([4])] depend on 9 and, hence, the GS energy E 
is also dependent on 6. 

To study the properties of the quantum paramagnetic 
phase existing in the vicinity of J2 = J\ /2 as well as the 
phase transitions to that phase we will consider gener- 
alized susceptibilities Xf that describe the response of 
the system to certain "field" operator ^F i 17 i 18 i 20 i 22 i 23 i 27 
To calculate such a susceptibility xf we add to Hamil- 
tonian (fTJ) a field term F = 8 O, where O is an operator 
that breaks some symmetry of H and the coefficient <5 de- 
termines the strength of the field. Using the CCM with 
either the Neel or the collinear reference state we calcu- 
late the energy per site E(8)/N = e(S) for H + F, i.e., for 
the Hamiltonian of Eq. ([T]) perturbed by the additional 
term SO. The susceptibility xf is then defined as 



XF = 



3 2 e(5) 



d5 2 



(7) 



5=0 



For the considered quantum spin model we have to 
use approximations in order to truncate the expan- 
sion of S and S. We use the well elaborated LSUBji 
schem o 64 i 66 ' 68 i 70 ' 73 in which in the correlation operators 
S and S one takes into account all multispin correlations 
over all distinct locales on the lattice defined by n or 
fewer contiguous sites. For instance, within the LSUB4 
approximation one includes multispin creation operators 
of one, two, three or four spins distributed on arbitrary 
clusters of four contiguous lattice sites. The number of 
these fundamental configurations can be reduced exploit- 
ing lattice symmetry and conservation laws. In the CCM- 
LSUB10 approximation we have finally 29605 (45825) 
fundamental configurations for the Neel (collinear) ref- 
erence state. Note, however, that for the calculation of 



the stiffness (the susceptibilities) the twisted reference 
state (the modified Hamiltonian H + F) is less symmet- 
ric, which leads to more fundamental configurations. As 
a result we are then limited to LSUB8 approximation. 

Since the LSUBn approximation becomes exact for 
n — > 00, it is useful to extrapolate the "raw" LSUBn 
data to n — > 00. Meanwhile there is a great deal of ex- 
perience how to extrapolate the GS energy e and the 
magnetic order parameter M. Most successful are the 
parameter fits of the form A(n) = Aq + Ai(l/n) Vl + 
^(l/ro) V2 where the fixed leading exponents v\ and v-i 
may be different for the different quantities to be ex- 
trapolated. For the GS energy per spin e(n) = ciq + 
ai(l/n) 2 + a 2 (l/n) 4 is a reasonable well-tested extrapo- 
lation ansatz^i^i^iii^iSiiZ^iS An appropriate extrap- 
olation rule for the magnetic order parameter for sys- 
tems showing a GS order-disorder transition i a 56 i 58 ' 59 i 75 
M(n) = b + bxil/n) 1 ' 2 + b 2 {l/n) 3 / 2 . For the spin stiff- 
ness the extrapolation p s (n) = cq + Ci(l/n) + C2(l/n) 2 
has been found to be reasonable^ Finally, for the sus- 
ceptibility we have tested several fitting functions, and 
we have found that the best extrapolation is obtained by 
the same fitting function xp{n) = co + ci(l/n) + C2(l/n) 2 
as for the stiffness. To check the reliability of this extrap- 
olation scheme we have also performed an extrapolation 
of the energy e(S) to n —> 00 by using the extrapolation 
formula e(n) = ao + ai(l/n) 2 + a 2 (l/n) 4 (see above) and 
a subsequent calculation of xf according to Eq. ([7]) using 
the extrapolated energy. We found that the deviations 
between both schemes are very small. 

In summary, the CCM approach automatically implies 
the thermodynamic limit N — > 00 (that is an obvious 
advantage in comparison with ED). However, we need to 
extrapolate to the n — ► 00 limit in the truncation index n, 
which is an internal parameter of the approach. Since no 
general theory is known how the physical quantities scale 
with n, we have to use extrapolation formulas based on 
empirical experience. Another feature of many approxi- 
mate techniques (but not of ED) is that they are based 
on reference states explicitly breaking some symmetry of 
the Hamiltonian. Although CCM also starts from a ref- 
erence state related to a particular magnetic LRO, it has 
been demonstrated that the CCM provides precise re- 
sults for the GS energy even in parameter regions where 
the magnetic LRO (i.e., the magnetic order parameter M 
calculated within CCM) vanishes ^^^^^SiH xhis is 
again an advantage of the CCM approach. 



III. GROUND-STATE PHASE DIAGRAM 

A. Ground state energy, magnetic order parameter 
and spin stiffness 

As already mentioned in Sec. |T1 the considered J\—Ji 
model has two semiclassical magnetic GS phases (small 
and large J 2 ) separated by nonmagnetic quantum phase 
(intermediate J2). To detect the quantum critical points 
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FIG. 2: (Color online) The GS energy per spin as function 
of J 2 obtained by CCM-LSUBn with n = 4, 6, 8, 10 and its 
extrapolated values to n — > oo using the extrapolation scheme 
e(n) = a + ai(l/n) 2 + a 2 (l/n) 4 . ED results for AT = 32 are 
shown by circles. 



by the above described CCM we discuss the magnetic 
order parameter M [see Eq. ([3], and the spin stiffness 
p s [see Eq. ©]. Both, M and p s , are finite in the mag- 
netically ordered phases but vanish in the intermediate 
quantum paramagnetic phase. 

For completeness, we show first the CCM and the ED 
GS energies per spin, e = E/N, in Fig. H The CCM 
curve consists of two parts corresponding to the Neel and 
collinear reference states, respectively. The dependence 
e{J%) for ED and CCM is qualitatively the same; how- 
ever, due to finite-size effects, the ED curve is below the 
CCM curves. Let us mention again that CCM GS en- 
ergy corresponding to the Neel (collinear) reference state 
is expected to be precise also in the intermediate quan- 
tum paramagnetic phase if J 2 is not too far beyond the 
transition points. 

Next we consider the magnetic order parameter in 
dependence on J 2 , sec Fig. [3] Note again that only 
for the magnetic order parameter M and the GS en- 
ergy we are able to solve the CCM-LSUBn equations 
up to n = 10, while for the stiffness and the suscepti- 
bilities we are restricted to n < 8. Hence the extrapo- 
lation to the limit n — > oo is most reliable for M and 
the estimation of the phase transition points by using 
the data for M is most accurate. The extrapolation 
to n — > oo shown in Fig. [3] is based the extrapolation 
scheme M(n) = b + b 1 (l/n) 1 / 2 + fe 2 (l/n) 3 / 2 and uses 
CCM-LSUBn data with n = 4,6,8,10. We find for the 
phase transition points between the semiclassical phases 
and the quantum paramagnetic phase Jf 1 = 0.447Ji 
and Jf 2 = 0.586Ji. To check the robustness of this ex- 
trapolation we have also extrapolated M using the data 
of n = 2,4,6,8,10 which leads to Jf 1 = 0.443 Ji and 
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FIG. 3: (Color online) Magnetic order parameter M versus 
J 2 obtained by CCM-LSUBn with n = 4,6, 8, 10 and its ex- 
trapolated values to n — > oo using the extrapolation scheme 
M(n) = bo + 6i(l/n) 1/2 + fe 2 (l/n) 3/2 . 



Jf 2 = 0.586Ji. Those values Jf 1 and Jf 2 are in agree- 
ment with CCM predictions of Rcfs. [56l andl58l 

Although the behavior of the extrapolated values of the 
magnetic order parameter around Jf 1 and Jf 2 presented 
in Fig. [3] shows a continuous behavior near Jf 1 an d near 
Jf 2 , it is obvious that the decay of the collinear order 
parameter to zero at Jf 2 is much steeper than the decay 
of the Neel order parameter at Jf 1 . That might give some 
hint of a first-order phase transition from the collinear 
to the paramagnetic phase, in contrast to a continuous 
transition from the Neel to the paramagnetic phase.— i 16 i 20 

In addition to the magnetic order parameter, another 
way to find the phase transition points is to consider the 
spin stiffness p s which is nonzero in a magnetically long- 
range ordered phase but vanishes in the magnetically dis- 
ordered quantum phase. The spin stiffness measures the 
distance of the ground state from criticality^i and consti- 
tutes together with the spin- wave velocity the fundamen- 
tal parameters that determines the low-energy dynamics 
of magnetic systemsi 78 ' 79 i 80 In order to calculate the stiff- 
ness directly using Eq. © we have to modify both the 
reference (Neel and collinear) states by introducing an 
appropriate twist 9; see Fig. Q] The CCM LSUBn re- 
sults for spin stiffness as well as the extrapolated values 
for both reference states as a function of J 2 are given in 
Fig. |U The results show that approaching the magnet- 
ically disordered phase the stiffness is decreased until it 
vanishes at J 2 = 0.466 Ji coming from the Neel phase and 
at J 2 = 0.578 Ji coming from the collinear phase. These 
values obtained by extrapolation including up to LSUB8 
data are in reasonable agreement with the critical points 
determined by extrapolating M. Note that our data for 
p a are also in good agreement with correspon ding results 
of the other methods (see Refs. HUG! and H^. Note 
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FIG. 4: (Color online) The spin stiffness p a versus J2 ob- 
tained by CCM-LSUBn with n = 4, 6, 8 and its extrap- 
olated values to n — > 00 using the extrapolation scheme 
p a {n) = co + ci(l/ra) + C2(l/n) 2 . 



instead of series-expansion data we use CCM and ED 
data. Interestingly, we will arrive at a different conclusion 
concerning the nature of the phase transition. 

The second type of arguments supporting the first- 
order phase transition from the Neel to the paramagnetic 
phase were based on series-expansion data for several sus- 
ceptibilities that test a possible valence-bond solid (VBS) 
order in the paramagnetic phase. In what follows we use 
the CCM and ED to compute four different susceptibili- 
ties Xj defined in Eq. ((7|) for the J1-J2 model. The corre- 
sponding perturbations (fields) Fj = S Oj, j = 1, . . . , 4, 
are given by 

F, = 5 J2(XX^+i^ (8) 
F2 = S (s hj s l+1 j - SijSij +1 ) , (9) 

f 3 = s J2xr +3 (stj^+w + siA+w) > ( 10 ) 



further that similarly as for M we observe also for p s 
that the curvature near the critical points is different at 
Jf 1 and at J| 2 that might be again a hint of the different 
nature of both transitions. 

To summarize, the CCM results for the GS energy, the 
magnetic order parameter, and the spin stiffness support 
a general physical picture known from earlier numeri- 



cal studies (including ED 



3.4.6.9.10 



variational quantum 



< J2 < 



Monte Carlo , 17 ' 21 and series expansion o 22 ' 27 ) . For in- 
termediate values of J2, 
(0.44 ± 0.01) Ji and Jf 
magnetic order. 



TCl 

J 2 



jc2 



with J| x 



(0.59 ± 0.01) Ji there is no 



B. Order of the phase transition: Generalized 
susceptibilities 

While the phase transition from the collincar to the 
paramagnetic phase is most likely of first order ) 6 ' 16 ' 20 
concerning the nature of phase transition from the Neel 
to the paramagnetic phase so far no conclusive answers 
are known. However, the question about the order of the 
phase transition from the Neel to the paramagnetic phase 
is of great interest in particular in connection with the 
validity of the Landau- Ginzburg paradigm. 40 ' 41 Very re- 
cently a number of arguments by Sirker et al.— based on 
series expansions and spin-wave theory were given that 
this transition is of first order. We reconsider this issue 
below using CCM and complementary ED results. 

The first type of arguments in favor of the first-order 
phase transition from the Neel to the paramagnetic phase 
presented in Ref . was based on the combination of field 
theory with series-expansion data. In what follows we 
use the same approach as that of Sirker et al.'^l however, 



(11) 



1 -J 



where i,j are components (integer numbers) of the lat- 
tice vectors of the square lattice [see Fig. [5]), where we 
visualize perturbation terms |8j - (fTT|) ]. The above def- 
initions, Eqs. ijHJ) - (ITT|) . arc in accordance with previ- 
ous discussion o 17 ' 20 ' 22 ' 23 and^ 7 - of possible valence-bond 
states or broken symmetries in the magnetically disor- 
dered quantu m p hase. Previous results for \i can be 
found in Refs. Hlili2(23| and 
22 and 



27, for X2 in Refs. [17| and 
27. Note that in Refs. HI 



H2I120I12: 
1271 . and for X3 in Refs 

and [27l the results for the perpendicular \3 [i- e -, the field 
F3 = SOz contains only x and y components, see Eq. (|10l) ] 
were reported only. For reasons of comparison with the 
available series-expansion data we consider in the present 
study also the perpendicular X3- To our best knowledge 
so far no data for the susceptibility %4 are published^ 2 - 
Note that all susceptibilities defined by Eqs. ([5]) - (jTTj) 
break the symm etry of the initial square lattice, for de- 
tails, see Refs. [THaHH and EH The susceptibilities xi 
and Xi are most interesting, since they belong to order- 
parameter operators Oi and O4 probing directly possible 
valence-bond ordering. As discussed in Ref. [4l| they can 
also be interpreted as a single complex order parameter 
with a different phase for the two patterns. Note that for 
the field F\ we have chosen the x-axis for the alignment 
of modified nearest-neighbor bonds, see Fig.[S^,. Alterna- 
tively, the y-axis can be chosen. It is worth mentioning 
that the field F4 [Eq. [IT] is a sum of fields F\ aligned 



and hence 



along x and y axes, i.e., F4 — ± x -r j- , 

(x) (y) 

Xi = x\ + Xi ■ If, in addition, a symmetry with re- 
spect to a 7r/2-rotation in the square-lattice plane holds 
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FIG. 5: (Color online) Illustration of perturbations (fields) 
Fj related to generalized susceptibilities \j '■ ( a ) perturbation 
Fi ©, (b) perturbation F 2 ©, (c) perturbation F 3 JTO]) and 
(d) perturbation Fi Dark (red) [light (green)] shadows 

correspond to enforced [weakened] exchange couplings. 



(that is, however, not the case, e.g., for the CCM calcu- 
lations for large J2), one has = Xi an d Xi — %Xi- 

Analyzing the behavior of the susceptibilities as J2 ap- 
proaches the critical value J^ 1 we will again arrive at a 
different conclusion in comparison to that in Ref . I27I . 

We begin with the examination of the order of the 
phase transition from the Nccl to the VBS state. In con- 
trast to the transition from the VBS to the collincar 
state where an energy level crossing indicates a first- 
order transition ) 16 ' 27 ' 56 the energy behaves smoothly as 
J2 varies around Jf 1 and a more sensitive method for 
distinguishing between first- and second-order transitions 
has to be applied^ For that we consider the GS energy 
e(S) for Hamiltonian ((T|) perturbed by the field F-y = SOi 
[Eq. We have performed CCM calculations for e(S) 
choosing the Neel state as the reference state and extrap- 
olating LSUBn data with n = 4, 6, 8 according to the 
scaling law e(n) = a + ai(l/n) 2 +a2(l/n) 4 [see Fig.[S^]. 
We have also performed complementary ED for a finite 
square lattice of N — 32 sites [see Fig. [BJd] for a qualita- 
tive check of the CCM results. The obtained dependence 
e(S) may be fitted for a fixed J2 to the following polyno- 
mial form 



a -2 , b xi , c. 6 



e(6) e(0) = + + -5 



(12) 



To determine the order of the phase transition we use 
the method described in Ref. H3- For a two-dimensional 
antifcrromagnet, the system could be described by the 
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FIG. 6: (Color online) The GS energy e(S) — e(0) versus 
square of field strength 5 for H + SOi [see Eq. J8], for 
J 2 = 0.0, 0.1, 0.2, 0.3, 0.4 and 0.45 (from top to bottom), 
(a): CCM results extrapolated to n — > 00 using the extrapola- 
tion scheme e(n) = ao + ai(l/n) 2 + a2(l/n) 4 . (b): ED results 
for N = 32. The displayed curves might be compared to the 
ones in Fig. 1 of Ref. [23 where corresponding series-expansion 
data for e(S) are reported (however, only up to J2 = 0.3). 



following 0(3) model: 



H, = i 
2 



{d t vf + c z v (vv)' z + m 2 v v 2 



4 V ; 



(13) 



Consider now the case that we are in the magnetically 
ordered phase and add the field F\ [Eq. [8] with \5\ <C 1. 
The Neel order will then coexist with a small dimcriza- 
tion described by a scalar field 



H 



4> 



(dt^ + c^r+m 2 ^ 2 



"T 



6 



(14) 

The fields v and 4> are not independent, and the interac- 
tion between them reads 



Hi, 



2 



(15) 
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FIG. 7: (Color online) The coefficient b of the quartic term 
in Eq. (frj]) obtained from a fit of the CCM data in Fig. [6k 
and the ED data in Fig. [6b in dependence on J2. This figure 
might be compared to Fig. 3 of Ref . I27I. Inset: the coefficient 
b versus J2 shown for small J2 with an enlarged scale. 

The effective field theory in the ordered phase for S 7^ is 
then given by H = H v + H^ + Hi nt . Combining Eqs. (fl"3]) - 
(I15[) we will have a nonzero GS expectation value 

W = ^-^ 3 + ^^ 5 + ^ 7 ), (16) 

with A = + u v< j)(v) 2 . Equation ([16]) leads to a GS 
energy given by 

<*>-<* = °) = -^ 52 + ^ 4 + ^^ 56+0 ^)- 

(17) 

The coefficient of the S 4 term in Eq. (fT7|) may be positive 
or negative depending on the sign of the parameter u<f,. 
In the case of > we have a second-order transition 
with respect to at a critical point, and a first-order 
transition if < 0. 

Using the polynomial in Eq. (fT2")) we have fitted the 
data of e(<5), <5 2 = . . . 0.09 for various J 2 including values 
near the critical point Jf 1 [see Fig. [6k,] ■ We find that the 
coefficient of the 5 4 -term b is negative for small values of 
J 2 but becomes positive if J 2 approaches Jf 1 ; see Fig. [7J 
This behavior is found for the CCM data as well as for 
the ED data. In particular, b calculated by the CCM 
(calculated by the ED) changes its sign at J2 « 0.35 (at 
J 2 «0.31). 

Comparing Fig. [7] with the results reported in Fig. 3 
of Ref. [27] we note that CCM data for J2 below 0.2 are 
in reasonable agreement with series expansions, linear 
spin-wave theory, or mean field spin-wave theory [in par- 
ticular, CCM yields b(J 2 = 0.1) « -0.40, b(J 2 = 0.2) « 
-0.86, b{J 2 = 0.25) w -1.29, b(J 2 = 0.3) « -1.73 that 
is in between the series-expansion data and the spin- 
wave theory results]. A drastic difference between the 



series-expansion data and the CCM results emerges if J 2 
approaches the critical value Jj : The series expansion 
gives b < whereas the CCM and the ED yield b > 
for J 2 — > J 2 l . We recall that any predictions from spin- 
wave theory for the considered J\-J 2 model are likely to 
be unreliable if J 2 exceeds 0.354 Combining Eqs. (fl"2"l) 
and Q17|) we get b = u^a 4 and determining a and b us- 
ing CCM data (Fig. EJi) for J 2 = 0.36 . . . 0.42 we find 
w w 0.75 > 0. 

In summary, the presented CCM and ED data, in con- 
trast to series-expansion data of Ref. I27I do not support 
a weak first-order phase transition from the Neel to the 
VBS stated! but give evidence that this transition is con- 
tinuous. 

Next we examine the susceptibilities associated with 
probing fields © - (HU) directly. The CCM results are 
shown in Fig. [HI We also present in this figure the ED 
data for N = 16, 24, and 32 lattice in the insets. (We 
do not show N = 24 results for xi and \2 since the sys- 
tem of rectangular shape perturbed by F\ or F 2 does not 
possess symmetry with respect to a 7r/2 rotation in the 
square-lattice plane.) Note that a sophisticated finite- 
size analysis has to be performed in order to derive the 
correct behavior of susceptibilities in the thermodynamic 
limit iil Such an analysis goes beyond the scope of the 
present study since we use the ED data as a qualitative 
check of our CCM results only. We notice here that al- 
though xi an d X4 are related to each other (see above), 
they are calculated completely independently. We have 
confirmed the expected relation between these suscepti- 
bilities thus providing an additional double check for our 
numerics. 

As it has been already mentioned above, the suscep- 
tibilities xi, X2j and X3 were calculated in earlier stud- 
ies using different methods. Our CCM results for xi 
and xi are m a good quantitative agreement with series- 
expansion results reported for J 2 = . . . 0.5 in Refs. I20I 
and I27I. [For instance, one can compare our CCM data, 
l/Xi(J 2 = 0.3) w 0.92, l/xi(J 2 = 0.35) w 0.66, and 
X2 (J 2 = 0.3) w 0.90, X2(J2 = 0.35) w 1.06, with the data 
shown in Figs. 2 and 5 of Ref. H3j The CCM results for 
Xi and %2 also qualitatively agree with variational quan- 
tum Monte Carlo method and ED results reported (for 
some J 2 only) in Ref. QjJ The CCM results for X3> how- 
ever, exhibit a different qualitative dependence on J 2 as 
J 2 approaches Jf 1 in comparison with series-expansion 
datai^ZZL Compare, e.g., Fig. 6 of Ref. and Fig. [St 
of the present paper. According to series-expansion data 
X3 decreases by about 20% as J 2 increases from to 0.4. 
In contrast, according to CCM data shown in Fig. [5J: X3 
increases by a factor of about 4 as J 2 increases from to 
0.4. 

Let us now discuss some general features of the gen- 
eralized susceptibilities shown in Fig. [5] Obviously, a 
divergence of a certain susceptibility (or 1/x — > 0) at 
a particular value of J 2 indicates an instability of a GS 
phase regarding to a possible different GS order. It can 
be seen from Fig.[5J that all susceptibilities increase with 
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FIG. 8: (Color online) The inverse susceptibilities (a) 1/xi, 
(b) l/x2, (c) 1/X3 (please note the scaling factor 0.01 at the 
j/-axis), and (d) 1/X4 versus J2 obtained within the CCM 
LSUBn approximation with n — 4, 6, 8 and extrapolated to 
n — > 00 using x{ n ) ~ c o + ci(l/n) + C2(l/n) 2 . Insets: the 
same as in the main panels but using ED for finite lattices of 
N = 16, 24, and 32. Panel (a) might be compared to Fig. 2 of 
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FIG. 9: (Color online) Susceptibilities 1/xi (red) and 1/X4 
(gray) around the critical point JJ 1 . Bold curves correspond 
to the CCM curves shown in Figs. |8^ and (8ji. Thin lines 
obtained from a linear fit of the CCM data for < J2 < J| • 
Extrapolated (thin) lines become zero at J2 w 0.47. 



growing J2 in the Neel phase. Near the critical point Jf 1 
both 1/xi and l/x4 (CCM data imply X4 = %Xi) are sig- 
nificantly smaller than I/X2 and 1 /\3 , indicating that the 
valence-bond states belonging to the columnar dimerized 
and plaquette patterns are favorable in the magnetically 
disordered quantum phase. 

A similar behavior of xi and Xi (CCM data imply 
X4 = xf^ + Xi) i s observed if J 2 approaches Jf 2 form 
the collinear phase, i.e., from J 2 > Jf 2 . On this side the 
behavior of X2 and X3 is not conclusive, since both are 
already large in the collinear phase. 

The behavior of the susceptibilities xi and X4 ( = 2xi) 
near the critical point Jf 1 is shown in more detail in 
Fig. [5J Obviously, approaching Jf 1 from the Neel phase, 
Xi ixt) becomes very large; it however, remains finite in 
the region around J^ 1 up to J2 = 0.55. That might be 
attributed to limited accuracy of CCM results since: (i) 
we have data only up to LSUB8 for extrapolation and 
(ii) and LSUBn data based on the Neel reference state 
may become less accurate for values of J2 exceeding Jf 1 . 
However, if the phase transition with respect to the corre- 
sponding VBS order parameter characterizing the quan- 
tum paramagnetic phase would be of second order we 
may expect an almost linear decrease in the inverse sus- 
ceptibility if J2 approaches Jf 1 , i.e., 1/xi oc (Jf 1 — J2) 7 * 
with 74, ss liSI Hence a linear fit of the CCM data of 
1/Xi (I/X4) versus J2 using data points only within the 
Neel ordered region < J2 < Jf 1 might give reason- 
able results. We find that the linear fit for 1/xr (I/X4) 
vanishes at the point J2 « 0.47Ji; see Fig. [9j This is in 
agreement with the scenario of deconfined criticality that 
predicts such divergence if the deconfined critical point 
is approached from the Neel phase. 27 i 41 
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To conclude this part, the CCM and ED data for all 
examined susceptibilities, xi> \2, X3i an( i Xii exhibit 
an enhancement while the system runs out of the Neel 
phase. This enhancement is most pronounced for xi ix^)- 
Moreover, \i (Xi) diverges at a value of J2 close to the 
quantum critical point Jf 1 « (0.44 ± 0.01) J\ determined 
by the most accurate data for the Neel order parameter 
M . This finding is consistent with the predictions for 
a deconfined quantum critical point 4^ Furthermore we 
find that our CCM data for xi an d X2 agree with the 
series-expansion data^i In contrast, for X3 we observe 
a qualitatively different behavior. Finally, the enhance- 
ment or divergence of the considered susceptibilities if J2 
approaches Jf 1 indicates that the translational symme- 
try is broken in the quantum paramagnetic phase, i.e. 
most likely a spatially homogeneous spin-liquid phase for 
Jf 1 < J2 < J£ 2 can be excluded. 

IV. SUMMARY 

To summarize, in this paper we have applied the CCM 
in high orders of approximation to the spin-1/2 J\-Ji 
Heisenberg antifcrromagnct on the square lattice and 
present a comprehensive analysis of the GS phase dia- 
gram of the model. For this purpose we have calculated 
the GS energy, the magnetic order parameter, the spin 
stiffness and several generalized susceptibilities. Our re- 
sults enrich the list of available data and are comple- 
mentary to other existing results obtained using differ- 
ent approximate methods such as series expansions or 
variational quantum Monte Carlo for the spin-1/2 J1-J2 
square-lattice Heisenberg antifcrromagnct. In addition 
to the CCM results we present also ED results that are 
found to be in good agreement with the CCM data. 

Our findings confirm the basic picture discussed earlier: 
For intermediate values of Jf 1 < J2 < Jf 2 the ground 



state is a paramagnetic quantum state. The CCM pre- 
diction for the boundaries of the paramagnetic region is 
Jf 1 w (0.44 ± 0.01) Ji and Jf 2 w (0.59 ± 0.01) Ji. To 
discuss the nature of the quantum phase transition from 
the semiclassical Neel phase to the quantum paramag- 
netic state at Jf 1 we use the CCM (and ED) data as an 
input for the method developed in Ref. H3 to distinguish 
between a first- and a second-order transition. Our anal- 
ysis leads to the conclusion that the phase transition from 
the Neel to the paramagnetic state at Jf 1 is second order. 
This outcome contradicts the conclusion of Ref. [27] based 
on series-expansion data, but agrees with the deconfined 
critical point scenario proposed in Refs. HO andl-Ltl, 

Another way to check the predictions of the theory of 
deconfined quantum criticality is to examine the suscep- 
tibilities related to order parameters of a possible VBS 
ordering emerging, if the critical point is approached from 
the magnetically ordered Neel phase. The obtained data 
shown in Figs. [8] and [9] give another hint that xi {Xi) 
diverges at Jf 1 which does not contradict the deconfined 
critical point scenario . 40 ' 41 

Finally, the divergence or enhancement of the general- 
ized susceptibilities obtained by CCM and ED approach- 
ing Jf 1 from the Neel phase gives evidence in favor of 
ground states breaking translational symmetry. There- 
fore, our data yield further arguments against a struc- 
tureless (i.e., a spatially homogeneous) spin-liquid state 
without any LRO. 
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